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Abstract. The NA61/SHINE collaboration performed measurements of pC interactions at 
31 GeV/c beam momentum with a full size replica of the T2K target (1.9 interaction length) 
during a pilot run in 2007. Larger statistics runs were also conducted in 2009 and 2010. The 
NA61/SHINE setup consists in a large acceptance spectrometer located on the H2 beamline of 
the SPS at CERN. For the first time, the kinematical phase space of interest for an accelerator 
based neutrino experiment (i.e. kinematical phase space of pions/kaons exiting the target and 
producing neutrinos in the direction of the near and far detectors) is fully covered by a single 
hadron production experiment. In a first stage, yields of positively charged pions were measured 
at the surface of the target. The analysis of the 2007 data set presented here demonstrates that 

a) high quality long target data were successfully taken with the NA61/SHINE apparatus, and 

b) for the first time, the T2K neutrino flux predictions can effectively be re-weighted with the 
NA61/SHINE long target data. 

Invited paper to NUFACT11, Xlllth International Workshop on Neutrino Factories, Super 
Beams and Beta Beams, 1-6 Aug. 2011, CERN and University of Geneva, submitted to IOP 
conference series. 

1. Motivations for long target measurements for T2K 

A prediction of the neutrino beam content at the T2K (pQ, [2]) far detector is given in Table 
[TJ The flux is predominantly (~ 95%) produced by the decay in flight of positively charged 
pions focused by the magnetic horns of the beamline. While ~ 30% of the v e flux is produced 
from the decay of positively charged kaons, ~ 50% of it is due to muons produced in the decay 
of the same charged pions that generate the v a flux. Thus, in a first stage, pion production data 
will constrain both f M and v e fluxes. 

Table 1. Composition of the neutrino beam and its various species at the far detector. 
Integrated values are quoted. Predictions obtained with the GCALOR model and horn currents 
set to 320 kA each (T2K was running with horn currents of 250 kA in 2010/11). 
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In terms of hadron production measurements, the pion contribution to the v„ flux at 
the far detector can be decomposed into a direct component and a target component. The 



first contribution originates from pions directly produced in the proton primary interaction 
(secondary pions) or in the decay of a secondary particle. The target contribution refers to all 
pions exiting the target or pions produced in the decay of an outgoing particle. The dependence 
of these contributions on the neutrino energy is depicted in Fig. [1] (left and middle panels). 
While the direct component contributes to ~ 60% of the flux at the beam peak energy 
of ~ 600 MeV, the target component contributes to ~ 90% of the flux. Thus, thin target 
measurements (p + C — > tt + + X at 30 GeV) for T2K constrain up to 60% of the flux 
prediction at the far detector, the remaining 40% being produced in secondary interactions in 
the target and elements of the beamline. Long target measurements would constrain up to 90% 
of the flux prediction. In this case, only 10% of the flux due to secondary interactions out of the 
target would require constraints from other data. 
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Figure 1. Ratio of (1-direct) [left] and (1-target) [middle] to total contribution for the flux 
at the far detector. Current error envelopes for the flux at the far detector [right] (for the 
analysis described in [4]). 



As depicted in Fig. [T] (right panel) , the total fractional error on the T2K flux prediction 
([3]) 01) & t the far detector is ~ 15% at the beam peak energy. Uncertainties arise from two 
main contributions: uncertainties related to the beamline setting (e.g. beam optics, beam 
direction, alignment and currents of the focusing horns) and those related to the re-weighting 
procedure based on the NA61/SHINE thin target data (within 5 to 10 %)[5j. The fractional 
error at the beam peak energy is dominated by uncertainties from the second contribution 
(i.e. pion/kaon multiplicity, secondary nucleon multiplicity and production cross sections). In 
particular, uncertainties on the re-weighting of secondary nucleon and production cross sections 
for secondary interactions are applied to ~ 40% (~ 10%) of the flux in the case of the thin 
(long) target based re-weighting. These uncertainties rely on interpolation between sparse data 
sets or comparison to models and are often poorly known. In the case where flux predictions are 
re-weighted by long target data, these contributions are taken into account at once and there is 
no need for error prone estimates of secondary interactions in the target. Long target data are 
therefore needed to obtain a precision of 5% or better on the absolute flux prediction. 



2. The NA61/SHINE long target measurements for T2K 

The 2007 NA61 /SHINE pilot data were taken with a replica of the T2K target at 30 GeV beam 
energy. The NA61/SHINE spectrometer consists in a set of 5 time projection chambers (TPCs) 
complemented by an array of time-of-flight (ToF) detectors located downstream of the TPCs, 
allowing for a full coverage of the forward pion production of interest for T2K (i.e. phase space of 
pions exiting the target and producing neutrinos in the direction of the near and far detectors). 
Details of the experimental setup are given elsewhere [5]. The trajectory of each beam proton 
track is reconstructed with a set of beam position detectors that allow to determine the position 
of the beam impact on the upstream face of the target. The usage of the ToF in the analysis 



assures that selected tracks in the TPCs originate from single proton interactions in the target. 
Data can therefore be normalized to the total number of protons on target (POT). 

Tracks reconstructed in the TPCs are extrapolated back through the magnetic field from their 
first measured point to the surface of the target. A point-of-closest-approach is found between 
the trajectory of the track and the surface of the target. This requires a precise knowledge of 
the relative alignment of the beam and target, which is determined by using both beam and 
TPC tracks. For the 2007 pilot run, the target was tilted with respect to the beam axis in 
both horizontal (~ 5 mrad) and vertical (~ 2.8 mrad) directions. The misalignment was taken 
into account in the extrapolation procedure and effects on the outgoing pion yields (additional 
systematic uncertainty) studied with dedicated Monte-Carlo simulations. The target is split into 
5 bins of 18 cm each along the beam (z axis), and a last bin corresponding to the dowstream face 
of the target that covers the very forward (below 40 mrad) pion production. Analysis cuts are 
optimized to improve the momentum and polar angle resolutions at the first measured point on 
track that determine the achievable resolutions on target, estimated to a z = 5 cm, a p /p = 3% 
and (Jo/ 9 = 6% for the longitudinal, momentum and polar angle resolutions respectively. As 
depicted in Fig. [2] (left panel), tracks can effectively be reconstructed on the surface of the target. 
The analysis coverage extends up to 20 GeV/c in the forward production region (< 40 mrad) 
and covers angles up to 280 mrad at lower momentum. 




Figure 2. Distribution of the point-of-closest-approach on target in the x-z projection [left], 
data distribution in the {ToF,dE/dx} phase space for 40 < 6 < 100 mrad and 2.4 < p < 3.2 
GeV/c [middle], and projection of the fit result along the ToF axis [right]. 

The particle identification (PID) relies on the energy loss (dE/dx) measurement in the TPCs 
and the time-of-flight that is used to compute the particle mass squared. The combination of the 
two signals provides a powerful PID over a large momentum range. Yields of positively charged 
pions exiting the target are extracted in bins of (p, 8, z) using a 2d-likelihood minimization (see 
Fig. [21 middle and right panels). Fig. [3] shows that such yields can effectively be measured at 
the surface of the target. 

The FLUKA simulation package [7j used in the T2K beam simulation was interfaced to the 
NA61 /SHINE simulation chain for which a PID based on data parametrizations was furthermore 
implemented. Data and Monte-Carlo are processed identically with the same reconstruction 
chain and PID analysis. Thus, reconstructed yields of positively charged pions can be compared 
to reconstructed Monte-Carlo predictions prior to any corrections necessary to get absolute yields 
(e.g. all Monte-Carlo based corrections such as acceptance and decay corrections, reconstruction 
efficiency, etc). Such a comparison is depicted in Fig. [3] (top left panel) with the corresponding 
re- weighting factors (top right panel), defined as the ratio of data to Monte-Carlo yields. Current 
uncertainties on the re-weighting factors are dominated by the statistical uncertainty which is 
typically 10% due to the very poor statistics of the 2007 pilot data. Systematic uncertainties are 



within 5-10%, the dominant source being PID (1-5%), followed by target misalignment (3%) and 
normalization (1.4%). Systematics from ToF efficiency, beam momentum and target density are 
all below 3%. 

Comparisons of FLUKA-standalone predictions for the T2K beam simulation to the 
NA61/SHINE long target data are also shown as an example in Fig. [3] (bottom panel), for 
predictions re-weighted with the NA61/SHINE thin target and long target data. More details 
on the NA61/SHINE long target data analysis and re- weighting procedure are given in [BJ. 
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Figure 3. Data and Monte-Carlo yields of positively charged pions normalized to the bin 
size and to the number of protons on target over the second longitudinal bin of the target for 
40 < 9 < 100 mrad [top left], and corresponding re- weighting factors [top right]. Comparisons 
of nominal FLUKA-standalone predictions with thin target and long target re-weighting to the 
NA61/SHINE long target data for different angular intervals [bottom panel]. 



Statistical uncertainties will be reduced with the 2009/2010 long target data sets, and further 
improvements are expected (in particular for the target alignment) to reduce the systematic 
uncertainties down to 5% or better. The pilot analysis of the 2007 NA61/SHINE long target 
data demonstrates that the T2K neutrino flux predictions can effectively be re-weighted with 
long target data. A first preliminary re-weighting of the Un flux prediction at the far detector 
was already implemented and compared to the prediction based on the 2007 NA61/SHINE thin 
target data ([5], [6]). 
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